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Abstract 



o 



A measurement of the decay in time of nuclei excited by an intense short laser pulse of energy Eq yields the Fourier 
transform of the autocorrelation function of the associated scattering matrix. We determine the optimal length (in 
time) of the pulse and evaluate the time-decay function using random-matrix theory. That function is shown to contain 
information not otherwise available. We approximate that function in a manner that is useful for the analysis of data. 
For Eq below the threshold energy E n of the first neutron channel, the time-decay function is exponential in time t while 
■it is the product of an exponential and a power in t for Eq > E n . The comparison of the measured decay functions in 
both energy domains yields an unambiguous and novel test of random-matrix theory in nuclei. 
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1. Purpose 

ELI, the "Extreme Light Infrastructure" , an ambitious 
European project to generate laser beams of extremely 
high intensity, is close to construction. Parallel to that 
development, experiments that will use those beams are 
being planned, and theorists are called upon to develop 
the concepts and tools needed for their analysis. At the 
workshop on ELI held in Palaiseau (France) April 27/28, 
2009, photonuclear reactions induced by an intense laser 
pulse of high energy (several MeV) received much atten- 
tion, see Ref. [if. The required high-energetic directed 
pulsed gamma rays are supposed to be produced by Comp- 
ton backscattering of a short laser pulse with much lower 
energy on a sheet of electrons ejected from a thin foil hit 
by the ELI pulse. 

In this paper, we present a theoretical study of nuclear 
reactions induced by short laser pulses of several MeV en- 
ergy. In particular, we address the following questions: (i) 
Which observable is measured in a photonuclear reaction 
induced by a very short laser pulse? (ii) Which novel in- 
formation is provided by data of that type? (iii) What is 
the optimal length in time of a short laser pulse for such 
a reaction? 

We focus attention on the main mode of nuclear excita- 
tion by gamma quanta of several MeV energy, the electric 
dipole mode. For a target nucleus in its ground state with 
spin J and parity it, dipole absorption leads to excited 
states with spin J ± 1 and opposite parity — 7r. The anal- 
ysis of experimental data will be simplest when the spins 
of the excited states are uniquely defined. Therefore, we 
consider an even-even target nucleus with a ground state 
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of spin zero and positive parity. Then the states excited 
by dipole absorption have spin 1 and negative parity. De- 
pending on excitation energy, these may lie below or above 
the first particle threshold. That is typically the thresh- 
old for neutron emission. In medium- weight and heavy 
nuclei, that threshold has an excitation energy of 5 to 8 
MeV. Nuclear states right above neutron threshold have 
been studied by time-of-flight spectroscopy in slow neu- 
tron scattering, mostly on even-even target nuclei. Such 
states appear as isolated s-wave resonances with spin 1/2 
and positive parity, with a typical spacing of 10 eV and 
a typical width of 1 eV. The statistical analysis of such 
resonances shows that the spacings and widths follow the 
predictions of random-matrix theory (RMT), see the re- 
view 0. In recent years, states in even-even nuclei with 
spin 1 and negative parity close to the first particle thresh- 
old forming the "pygmy dipole resonance" have been in- 
tensely studied experimentally, mainly with the help of the 
resonance fluorescence technique, see Refs. [3,IH and refer- 
ences therein. In Ref. Q deviations from RMT predictions 
were found that have so far not been fully explained the- 
oretically. Our work is based on the assumption that the 
states excited by an intense short laser pulse are governed 
by RMT. We explore the consequences of that assumption 
and propose an experimental test for it. 

For reasons given below, we consider a laser pulse of 
10 -19 to 10~ 20 s duration and a mean energy of several 
MeV. The pulse coherently excites a band of l - states. 
The band width is several 10 keV, the number of 1~ states 
involved is typically 10 3 to 10 4 . We refer to these states 
as to compound-nucleus (CN) resonances. These will sub- 
sequently decay. CN resonances below neutron threshold 
decay by gamma emission, those above neutron threshold 
preferentially by neutron emission. Gamma emission is 
possible but less likely. In both cases the detection of the 
emitted particle with highest energy unambiguously iden- 
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tifies decay into the ground state of the residual nucleus. 
We do not focus attention solely on that decay mode since 
in the case of gamma decay, the intensity may be too low, 
and a summation of intensities corresponding to several or 
many gamma decay modes leaving the target nucleus in 
its ground or one of its excited states, may be called for. 
In any case we deal with a two-channel situation. The 
incident channel 70 is defined by the target nucleus in its 
ground state plus a dipole gamma quantum, the exit chan- 
nel b by the emitted particle (b = n for neutron emission 
leaving the residual nucleus in its ground state, b — 7, 
with i = 0, 1, . . . , A for gamma emission leaving the resid- 
ual nucleus in its ground or any of A excited states). 

Formation and decay of the CN resonances being coher- 
ent processes, it is not possible to identify any particular 
CN resonance as the source of the emitted particle. In- 
terest in the data rather focuses on the time dependence 
of the decay. We answer the questions raised in the sec- 
ond paragraph by identifying the relevant observable and 
defining the optimal length in time of the laser pulse. We 
also work out the expected form of the time-decay func- 
tion. We do so by using results of the analytical approach 
to CN reactions developed in Ref. Q. 

Coherent resonance formation was previously addressed, 
see Ref. @ and references therein, the decay in time of 
CN resonances in Refs. @, H, 0. The deviations from 
the exponential decay in time predicted in Refs. 0, 0, 
Q have been observed in microwave billiards by Fourier- 
transforming the measured elements of the scattering ma- 
trix, see Refs. 13, U|, and the comprehensive summary 
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12j. Here we propose a direct measurement of the 
time-decay function in nuclei. We go beyond Refs. [zlH,® 
in addressing specifically the case of photonuclear reac- 
tions induced by an intense laser pulse. In Ref. [13[ the au- 
tocorrelation function of the total photodissociation cross 
section for a chaotic atom or molecule was studied. With 
the help of the optical theorem, that function is related 
to a two-point function similar in form to the expressions 
studied below. The approach was extended and general- 
ized in Refs. 
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2. Observable 

To describe the effect of a short light pulse, we con- 
sider the scattering wave function \&+ (E) . That function 
describes all reactions caused by a gamma quantum of en- 
ergy E with wave number k(E) incident on the target nu- 
cleus. Asymptotically, \P+ (E) has an incident wave with 
unit flux in channel 70 and outgoing waves in all channels. 
The amplitude of the outgoing flux in channel b is given 
by the element S b ^ (E) of the scattering matrix, and the 
cross section a bj0 (E) feeding final channel b is given by 



a b~f {E) - — \S b7o (E) 



(1) 



/ 



dEg(E)exp{ 



-iEt/h}*+(E) 



of scattering wave functions with different energies where t 
denotes the time. The envelope function g(E) is a smooth 
function of energy E centered at energy Eq (the mean 
energy of the laser pulse), has band width AE, and |g(£')| 2 
is normalized to unity. The amplitude of the outgoing flux 
in channel b is J AE g(E) exp{—iEt/h}S bl0 (E). The total 
scattered flux I bl0 (t) in channel b versus time is given by 



h l0 (*) 



dE 



9(E) 
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exp{-iEt/h}[S bl0 (E) - 6, 



&70. 



(2) 

Indeed, integrating the total scattered flux over time and 
using Eq. ((2j we obtain the energy-averaged cross section, 
/ dt/h I bl0 {t) = J dE \g(E)\ 2 a bl0 (E). In other words, 
Ib~f (t) gives the decomposition in time of the cross section 
induced by a laser pulse of band width AE and time length 
At = h/ AE. A short laser pulse is defined to have a band 
width AE that is large compared to the average spacing d 
of the CN resonances. We focus attention on times t that 
are large compared to At and on the associated long-time 
behavior 4°" g (i) of I bla (t). We do so because it may be dif- 
ficult to separate experimentally short-time contributions 
from the original laser signal, and we accordingly expect 
that the first round of experiments will focus on the long- 
time aspects of the reaction. The terms in Eq. @ that 
involve § bl0 are proportional to the Fourier transform of 
g{E), have time length At, and are, therefore, neglegted. 
Thus, 
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(3) 



1.2 



J dE g{E) exp{-iEt/h}S b ^ (E) 

^^JdE 1 JdE 2 g(E 1 )g*(E 2 )exp{i(E 2 -E 1 )t/h} 

x Sb 7o (Ei)S b7o (E 2 ) . 

Here and in what follows we use that k(E) changes lit- 
tle over the energy interval AE and replace k(E) by fco = 
k(Eo). (We note that J lon s still contains a short-time com- 
ponent, see Eq. (TI4")) below. That component will also be 
suppressed eventually). We relate ^°" s (i) to the S'-matrix 
autocorrelation function C bl0 (e) defined by 

C bl0 {e) = J dE \f(E)\ 2 S bl0 (E -s/2)SZ Jo (E + e/2) . (4) 

The average is performed with a smooth function \f(E)\ 2 
normalized to unity the range of which extends to infinity. 
An example is a normalized Lorentzian with a width that 
is very large compared to d. We deal with isolated CN 
resonances for which the average total resonance width T 
obeys T < d. Then the function C{, 70 (e) decreases rapidly 
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towards zero for |e| > d. In Eq. (|3]) we write E\ = E — e/2, 
E-2 = E + e/2. We use that g(E) is smooth over the 
correlation width of Cb 7o (e). Then g(E ± e/2) re g(i?), 
and Eq. ([3]) becomes 



(5) 
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The approximation leading to Eq. ([5]) requires AE J> d. 
Otherwise, finite-range-of-data errors have to be taken 
into account. If that condition is met, the long-time part 
of the signal if, 7o (t) measures the Fourier transform of the 
S-matrix autocorrelation function Cb 7o (e). The inequal- 
ity AE d yields the first constraint on the length in 
time of the laser pulse: At must be small compared to the 
Heisenberg time H/d. 

A second constraint arises because dipole absorption in 
nuclei is governed by the giant dipole resonance. Depend- 
ing on mass number the resonance occurs at excitation 
energies between 10 and 15 MeV and has a width T dip of 
several MeV. The resonance causes a secular variation of 
the average dipole strength of the CN resonances. For a 
clean theoretical analysis, it is desirable to separate that 
secular variation from the average over resonances taken in 
the first of Eqs. ©. That condition is met if AE < T dip . 

Combining the two constraints we find d >C AE <C 
r di P. With d re 10 eV and T di P re several MeV, the con- 
straints are met for AE in the range 10 keV to 100 keV, de- 
pending on excitation energy and mass number. That cor- 
responds to a time length At re 0.5 x 10~ 19 s to 0.5 x 10~ 20 
s. Under these conditions, irradiation of a target with an 
intense laser pulse yields information on the decay in time 
of the CN resonances. As shown in Refs. [1, that decay 
is not expected to be exponential in general. It would cer- 
tainly be exciting to measure directly the time-dependence 
of the decay of the CN, i.e., of the signal Ib l0 (t), see Sec- 
tion |31 even though measurements in microwave billiards 
have already provided clear evidence for non-exponential 
decay in these systems 1^, U, 12 1- Such measurements 



are possible because the average width for neutron decay 
is less than 1 eV and that for gamma decay even smaller. 
The overall decay time of the CN resonances is, thus, much 
longer than the duration of the incident laser pulse, and 
the two signals are clearly separated. 

We note that a band width AE re 10 ... 100 keV also 
guarantees that in exciting CN resonances right above neu- 
tron threshold, one avoids excitation energies where neu- 
tron emission leaving the residual nucleus in an excited 
state becomes possible. (The excitation energy of the first 
excited state in the residual nucleus is typically 100 to 
200 keV). This is desirable as otherwise in the expression 
for Cb~f (e) the additional neutron channel must be taken 



into account, even if only decay into the ground state is 
experimentally measured. 

The decay in time of CN resonances excited by a short 
and intense laser pulse is described by the Fourier trans- 
form of the S-matrix autocorrelation function Cb 7o (e). 
Is that information novel, or are other data available or 
within experimental reach that would yield the same in- 
formation? An obvious possibility would be, for instance, a 
measurement of the (70, vi) reaction cross section a ritl0 (E) 
from which the cross-section autocorrelation function 



J dE \f(E)\ 2 a nno (E + e/2)a nnn (E - e/2) . 



with f(E) as defined in Eq. ((4]) can be obtained. For iso- 
lated resonances, C" (e) is neither theoretically accessi- 
ble (in contrast to Cb 7o (e), see Section [3]), nor is it simply 
related to Cb 7o (e), see Ref. (l6| . Generally speaking, the 
coherent decay in time of the CN resonances is described 
by an amplitude correlation function. An intensity au- 
tocorrelation function such as C£ 7o (e) does not provide 
equivalent information (except in the Ericson regime where 
r ^> d). Thus, the information available from Cb 7o (e) is 
unique. 

3. S Matrix Autocorrelation Function 

The autocorrelation function for CN scattering has been 
calculated analytically for CN resonances that obey RMT 
statistics [5j. In Ref. [17] and references therein it was 
shown that this function correctly describes chaotic scat- 
tering. Without going into details, we summarize some 
salient features of that approach, see Refs. [Ill [13], and 
cite the result. The S-matrix with elements Sb a (E) is de- 
composed into an average part (Sba) and a fluctuating part 
Sf a (E), 

Sba(E) = (Sba) + Sl(E) . (6) 

The average is taken over the ensemble of random ma- 
trices. In what follows we use the equality of ensemble 
average and energy average which holds provided the lat- 
ter is performed over an energy interval containing a large 
number of CN resonances with a smooth averaging func- 
tion \f(E)\ 2 as in Eq. (g|. The average S'-matrix (Sba) 
describes the fast part of the reaction. Simple models like 
the optical model for elastic scattering or direct-reaction 
models yield reliable theoretical results for (Sba)- The fluc- 
tuating part S^ a (E) describes the slow part of the reaction, 
i.e., formation and subsequent decay of the CN resonances. 
Because of the complexity of the latter, the precise energy 
dependence of S^ a (E) cannot be predicted theoretically. 
Using a random-matrix model for the CN resonances one 
can, however, calculate the energy-autocorrelation func- 
tion of Sf (E). That function is given in terms of the av- 
erage spacing d of the CN resonances and of the elements 
of the average S'-matrix (S) which serve as input parame- 
ters. We first assume that (S) is diagonal and later discuss 
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modifications due to the non-zero non-diagonal elements. 
We have @ 

<SL(£i)(SL(£ 2 ))*> = (St(E-e/2)(St(E + e/2)y) 

2 p-\-oc pi 

= 11 / dA, / dA J 6a (A l7 A 2 ,A) 



I I Z7T£ 

x -/i(Ai,A 2 ,A) exp<{ -— (Ai + A 2 + 2A) 



(1 - T e A) 



(l + T e A 1 ) 1 /2(i +Te A 2 ) 1 /2 



(7) 



We note that the autocorrelation function © depends 
only on the energy difference e = E% — E\. The factor 
A*(Ai, A 2 , A) is an integration measure and is given by 



/i(Ai,A 2 ,A) 



(1 - A)A|Ai - A 2 | 



Ili=i[((l + A l )A l ) 1/2 (A + \i) 2 } ' 



(8) 



while 



Jfca(Ai,A 2 ,A) = (l + S ab )T a T b (9) 
Ai(l + Ai) 2A(1-A) 
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^(l+TaAiXl + T^Aj) (l-T a A)(l-T 6 A) 



<UT a 2 (l-T Q ) 
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A, 
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describes the dependence of the correlation function on en- 
trance and exit channels a and b. The product in Eq. © 
extends over all open channels e. The correlation func- 
tion © depends on the average level spacing d of the CN 
resonances and on the "transmission coefficients" 



T a = I - \(S aa }\< 



(10) 



Given (S aa ) and d, the right-hand side of Eq. (|Tj) is com- 
pletely known. The full S'-matrix autocorrelation function 
reads then as 



C ba (e) = S ab \(S aa )\ 2 + (SUe^SUEz))*} . 



(11) 



We turn to the case where (S) is not diagonal. This 
is of practical interest because the dipole operator gives 
rise to direct (7,, n) reactions and, thus, to non-vanishing 
elements (S nii ). We assume that (S n7i ) and (S nn ) are 
known and that (S li - yj ) is diagonal. (Non-diagonal con- 
tributions (5* 7i7j ) with i ^ j would be of second order 
in the electromagnetic interaction and, thus, negligible). 
We follow the work of Ref. [l8| summarized in Ref. 
If (S a b) is not diagonal, one has to determine the unitary 
transformation U ab that diagonalizes the transmission ma- 
trix P ab = S ab - Ec(Sac)(St c ) so that (UPU\ ab = S abPa . 
The 5-matrix transforms according to S — > S = USU T 
where (S) is also diagonal. The correlation function of S is 
given by Eqs. © to ©, with all T a replaced by p a . Since 
(S n ~ti) is governed by the electromagnetic interaction, we 



have that \(S nii )\ <C KSVm)! for all i. This allows us to 
use first-order perturbation theory to calculate U and the 
eigenvalues p a . We find that to lowest non- vanishing or- 
der in the electromagnetic interaction, the autocorrelation 
function of S = WSU* is equal to that of S as given in 
Eqs. © to ©, except for the replacement 



^ P"fi T^n I (Sfl'Ji ) I 

The autocorrelation function is then 



(12) 



C bl0 (e) = \(S b70 )\ 2 + (S^E^S^iEz))*) . (13) 

The term (S^ (Et)(S^ o (E 2 ))*) is given by Eqs. © to © 
with the replacement Eq. (fT2"j). 

We apply these results to the case of interest by choos- 
ing a = 70 and calculating the Fourier transform of C bla (e) 
and, from there, the intensity -^™ s (i)- We find (see Eq. ©) 

M/^( t ) = 27 rM(i)5 b70 K5 7070 )| 2 (14) 
de e^{iet/n}(S^ (E - e/2)(S^ (E + e/2))*) . 

The delta function in the first term on the right-hand side 
of Eq. ([Till signals that the contribution from (S) is instan- 
tanteous in time. That would hold for an infinitely short 
laser pulse. For the actual laser pulse, the signal will have 
the same time duration At as the pulse itself. In any case, 
that term does not contribute to the long-term behavior. 
Thus, 



^?(t) = c^ (t) 



(15) 



= J de e^et/h}(Sl {E~e/2){Sl {E + e/2)y). 

The first equation defines Cf, 7o (f). From Eq. © we have 
C bj0 (t) (16) 

= 2dT[ / dAi / dA 8{[td/(irh)] - [Ai + A 2 + 2A]) 
l=1 Jo Jo 

V I M(Ai, A 2 , A) g (1 + T J]i/^i + TiM i/2 



i=0 

/r, r I 1 A (1 ~ ^raA) . 

11 "' opcnJ + "' opcn (1 + T„A 1 ) 1 /2(i + T n A 2 )V2 > 
x Jb 7o (Ai, A 2 , A) . 

where the replacement Eq. (|12l) has to be made. The func- 
tion C{, 70 (t) gives the decay intensity of the CN resonances 
and is the object of central interest. The delta function un- 
der the integral in Eq. (|16[) shows that contributions due 
to the decay of the CN resonances are delayed (all three 
integration variables are positive). The symbol <5„. opcn is 
zero (unity) if the neutron channel is closed (open), re- 
spectively. 
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4. Decay in Time of the CN Resonances 

In this and the next Section we work out the time de- 
pendence of the time-decay function C& 70 defined in Eqs. (fT5 
and (|16p . with the replacement (|12[) . in a time domain 
where the signal is sufficiently strong for detection. The 
average correlation width T of the CN resonances is ap- 
proximately given by the Weisskopf estimate [l9j | . V = 
{d/(2Tr))J2 e p e . Since the CN resonances are isolated, we 
have Yld=oPn < The number (A + 1) of open gamma 
channels is very large, A ^> 1, so that p li <C 1 individually 
for all i. The value of the transmission coefficient p l0 can 
be obtained from the average total cross section for dipole 
absorption given by 



2-jr 



(17) 



According to the statistical model, (S' 7o7o ) is real and, for 
weak coupling to the channels, positive. From Eq. (1101) 
(with T replaced by p) we have {S Jo7o ) = v/l — P 7o an d, 
for p 7o <C 1, (5 7o7o ) « 1 — p 7o /2. That yields 



E( cr b7o) 



(18) 



For the transmission coefficient T„, we observe that the 
neutron has angular momentum one (zero) if the parity 
of the residual nucleus is positive (negative), respectively. 
We use the fact that for the CN resonances seen in slow 
s-wave neutron scattering, the average width T is about 
1 eV. Moreover, T is dominated by the neutron channel, 
T « T„. With T n = 2TiT n /d and d 10 eV that gives 
T n ps 0.6 for s-wave neutrons. The transmission coeffi- 
cient for p-wave neutrons is smaller by the p-wave angu- 
lar momentum barrier penetration factor kR (with R the 
nuclear radius and k the wave number). For a laser pulse 
of several 10 keV band width and a mean energy of 50 keV 
above neutron threshold, we have kR » 0.25, leading to 
T n rs 0.1 or 0.2. At an energy right above neutron thresh- 
old T n is considerably smaller, of course. Thus, T n is much 
larger than any of the p~ ti and of the same order as or even 
larger than Y^iPti- This shows that we must treat T n and 
p li in Eq. (|T6l) differently 

The evaluation of Eq. (|16[) seems to require the knowl- 
edge of all individual transmission coefficients p li for all 
photon channels. These are not known. However, Eq. (1161) 
can be much simplified so that it depends only on the to- 
tal decay width for gamma decay and on the transmission 
coefficients in the entrance, in the exit, and in the neutron 
channel. 

We are guided by the following observation. In Ref. [9j 
it was shown (see Eq. (6.12) of that reference) that if 

Tj <C 1, the Fourier transform of the S'-matrix au- 
tocorrelation function (S a b{E — e/2)S* b {E + e/2)) is given 

by 
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If the number of channels is large, that expression is ap- 
proximately given by exp[—(dt/h)(J2 i T i + 2T a + 2T b )] ) and 
the approximation is excellent for times t <C h/{dTi) for 
all i. If all Ti are approximately equal, that condition is 
met for all times for which the signal is detectable. That 
shows that for A> 1, the non-exponential decay predicted 
in Ref. |9| actually becomes unobservable: Deviations from 
the exponential decay form occur only for times for which 
the signal is too small for detection. We use that fact to 
simplify Eq. Q16p . We assume that the p li all have similar 
values and use the approximation (1— p li \)[(\+p li Ai)(l + 
p ri \ 2 )]~ 1/2 ~ exp(-[p 7i /2][Ai + A 2 + 2A]). That implies 



(l-P 7l A) 



i=0 



(I+^AO^I+Pt^) 1 / 2 



exp{— ~^^p 7i td/h} . 



(19) 

The arrow indicates that we have used the delta func- 
tion in Eq. ([TI)]) . We reiterate that the approximation 
Eq. (fTTJ]) is expected to be excellent for times that allow a 
detection of the signal although it may fail asymptotically 
(t — > oo). To interpret the right-hand side of expression 
Eq. (UJ]) we use that for T 7i « 1 we have T 7i = 27rr 7i /d 
so that {d/h)^2 i T~ ji = T^/h where T 7i and T 7 are the 
average partial and total widths for gamma decay of the 
CN resonances, respectively. That is the expected result. 
The replacement Eq. (|T2|) implies that the average partial 
widths and the total width for gamma decay are reduced 
by the direct reaction. This is plausible because some de- 
cay strength is taken away by that reaction. We denote 
the ensuing average total width for gamma decay by T 7 . 
The corresponding increase in neutron decay width is neg- 
ligibly small in comparison with T n . 

With the approximation leading to expression (|19[) . 
Eq. (TTBf becomes 

J As e^{ietlh}{Sl Q {E~e/2){Sl Q {E + e/2)y) 
= 2dcxp{-f 1 t/h}F bl0 (t) (20) 
where 

Fb 70 (t) (21) 



2 z'+oo pi 

dA 4 / dA 6([td/(irh)] - [A a + A 2 + 2A]) 
Jo 



=n 

x - / Lt(A 1 ,A 2 ,A) 



X ([1 ^n.opcn] "f" ^n,ope 
X Jfc 7[) (Ai, A2, A) . 



(1 - T n \) 



(i + r n A 1 ) 1 /2(i + r n A2) 1 / 



1/2 ) 



1 



The time-decay function in Eq. (|20l) is the product of two 
factors. The first is an exponential and describes the de- 
cay due to all gamma transitions (except for additional 
contributions from the entrance and the exit channels). 
The second factor Fb l0 (t), given in Eq. (f2"Tj). depends on 
the channels under consideration and on a small number 



of parameters: Via Eq. (|9]) F(, 7o (t) depends on the trans- 
mission coefficients in the entrance and exit channels and 
via the delta function on the average level spacing of the 
CN resonances. Moreover, Fb l0 (t) depends on whether the 
neutron channel is open or not, and on the value of T n . In- 
terest focuses on Fb l0 (t) because it gives rise to observable 
modifications of the exponential decay. 

We expect Fb l0 (t) to vanish for t < (this is confirmed 
by the delta function in Eq. (|2I[1 ). to rise to a maximum 
at some positive value of t, and to decay towards zero for 
t — > oo. Obvious questions are: At which value of t does 
the maximum of i 7 i, 7o (t) occur? How steep is the rise for 
small positive values of t? What is the form of the de- 
cay for values of t beyond the maximum? Some of these 
questions can be answered analytically. For small positive 
times, F\, l0 (t) rises quadratically. Indeed, because of the 
delta function in the integrand of Eq. (|2T]l both i 7 6 7o (i) 
and its first derivative vanish at t = 0. For larger values of 
t the behavior of Ff, 7o (t) is expected to differ for the three 
possible cases: (i) the neutron channel is closed, (ii) the 
neutron channel is open and neutron decay is measured 
(b = n) and (iii) the neutron channel is open but gamma 
decay is measured (b — ji). The results of Ref. [9J sug- 
gest that in case (i) the peak of F lil0 {t) is followed by a 
decay of the form t~ 2 . However, since both p l0 and p 7i 
are very small in comparison with T 7 , that decay is very 
slow, and the behavior of the time-decay function (f!?U|) be- 
yond the peak of F lil0 (t) is governed by the first factor in 
Eq. ([201) . i.e., is purely exponential. In case (ii) we expect 
that T n d/H is at least as large as T 7 . The large flux into 
the neutron channel should shift the peak of F n7o (t) to- 
wards smaller values of t than in case (i). The decay in 
time of F nia (t) beyond its peak should be governed by the 
neutron channel, too, and should asymptotically be pro- 
portional to t~ 3 / 2 . The decay time is comparable with or 
smaller than h/T 7 , and modifications of the exponential 
decay form should be detectable. Similar statements ap- 
ply in case (iii) except that now beyond its peak the time 
decay of F lil0 {t) is asymptotically given by rj -1 / 2 . 



5. Numerical Results 

Further insight into the time dependence of the time- 
decay function (|16|) is obtained by numerical simulation. 
Taking A = 49, choosing all p li equal to T, and using 
different sets of values for T and T n , we have found that 
Eq. (|20[) is in all cases an excellent approximation to Eq. (|I6D 
for those values of t for which the signal is detectable, both 
when the neutron channel is closed and when it is open. 
By way of example this is shown in Figs.[TJto|H with the 
black lines showing the analytic function Eq. (fT6"f and the 
red crosses the approximation Eq. (|20[) . That is an impor- 
tant result as it enables the analysis of data in terms of 
a few parameters and without knowledge of the individ- 
ual values of the transmission coefficients in every gamma 
channel. 
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Figure 1: Intensity C 7070 (t) versus time t in units of h/d. The neu- 
tron channel is closed, and A = 49 inelastic gamma channels are 
open, with p~ li = p 70 = T, i = 1, ...,49 and values of the transmis- 
sion coefficients chosen as indicated in the insets. Solid line (color 
online: black): Eq. 1 161 1 . Crosses (color online: red): Approximation 
Eq. J20I) . Open circles (color online: blue): Fit of an exponential 
cxp(— a,2t) to the data, resulting in f» JZj—o P-y, ■ 



The rise in time of the time-decay function C& 70 (t) is 
very steep in all cases. The function reaches its maxi- 
mum at times of order h/d. The maximum value of Ct> 70 
is mainly determined by the factor T a Tb in Eq. ([9]); that 
explains the enormous difference in scale in the figures. 
With increasing values of the transmission coefficients the 
maximum is shifted towards smaller values of t. These are 
very short times; it is not clear whether in the first round 
of experiments the decay signal can be clearly separated 
from the signal due to the short pulse itself (delta function 
in Eq. ([T4"])). Therefore, we have focussed attention on the 
decay in time of C& 7o (t) beyond its maximum. 
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Figure 2: Same as Fig.f^but for an inelastic gamma channel, C 7i70 (t) 
with i ^ 0. 

When the neutron channel is closed, C& 7o (i) is very 
well approximated by an exponential, as expected. This is 
shown in Figs. Q] and [5] with the exponential fit shown as 
blue open circles but applies equally to all other cases cal- 
culated. The exponential is the same for all gamma chan- 
nels. This is important because the signal is expected to 
be very weak for every single gamma channel. Summation 
over many such channels does not affect the form of the ex- 
ponential. Measurements of a restricted sum J2i C 7i7o (t) 
would yield the average gamma width of the CN reso- 
nances. 

Taken by itself, exponential decay is expected and not 
very exciting. The situation changes when the neutron 
channel is open. This is shown in Figs. [3] and SJ Because 
of the comparatively large value of T n the neutron yield 
is significantly larger than the yield in any single gamma 
channel. Moreover, the time decay function differs signifi- 
cantly from an exponential, both in the neutron and in the 
gamma channels. As indicated in the captions, the curves 
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Figure 3: Logarithmic plot of Ci, 7 „ (t) versus t (in units of h/d) for 
an open neutron channel and A = 49 inelastic gamma channels, with 
Pn = Pia = T, i = 1, ...,49, and values of the transmission coeffi- 
cients T = 0.0004, T n = 0.2 and of the final channel b as indicated in 
the insets. Solid line (color online: black): Eq. II16H for large times 
t. Crosses (color online: red): Approximation Eq. (12011 . Open circles 
(color online: blue): Fit of t ai exp(— <i2i) to data. In all three cases 
a 2 ~ 0.02 which is fa 50 x T. For b = n we find ai ~ 1.46 fa 3/2, for 
b = 7o> li we have a\ ~ 0.48 fa 1/2. 

were fitted with a function of the form t ai exp(— a^t) (blue 
open circles). The best-fit value of ai is approximately 
equal to the sum of the transmission coefficients of the 
gamma channels. The exponent a± agrees approximately 
with the result given at the end of the last Section. Both 
results are in agreement with our expectations. 
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Figure 4: Same as Fig. [3] but for a different set of transmission co- 
efficients T„ = 0.4, T = 0.0016. In all three cases a 2 ~ 0.077 which 
is fa 50 X T. For b = n we have ai ~ 1.39 fa 3/2, for b = 70, 7; we 
obtain a\ ~ 0.51 fa 1/2. 



6. Summary and Conclusions 

In nuclear reactions induced by short laser pulses of 
several MeV energy, the observable of interest is the time- 
decay function of the CN resonances. Provided the length 
At of the laser pulse is chosen optimally, the time-decay 
function is given by the Fourier transform of the S'-matrix 
autocorrelation function. For that to be true, At must be 
large compared to the Heisenberg time h/d and small com- 
pared to the width of the giant dipole resonance. This fixes 
At to values between 0.5 x 10~ 19 s and 0.5 x 10~ 20 s, de- 
pending on mass number. The time-decay function com- 
prises information on amplitude correlations of CN reso- 
nances which cannot be obtained from other observables. 

We have calculated the time-decay function under the 
assumption that the laser-excited CN resonances are de- 
scribed by random-matrix theory. Our Eq. (|20p gives an 
excellent approximation to that function. It depends on 
a small number of parameters only and is useful for the 
analysis of data. We have shown how to estimate these 



parameters from existing data (transmission coefficient for 
neutrons, average cross section for dipole absorption). The 
time-decay function rises steeply with time and reaches a 
maximum a short time after the initial laser pulse has hit 
the target. That time is of the order of the Heisenberg 
time. The further development in time of the time-decay 
function depends on whether the neutron channel is closed 
or open. In the first case, the time-decay function de- 
creases exponentially. The decay width is given by the 
average total gamma decay width of the CN resonances 
and can be determined from data. In the second case, 
the time-decay function can be fit by the product of an 
exponential (again determined by the average total width 
for gamma decay) and a power law. The exponent of the 
latter depends on whether the final channel is the neutron 
channel or a gamma channel. 

An experimental confirmation of our predictions would 
establish a new unambiguous test of random-matrix the- 
ory in nuclei. In addition, it would make it possible to 
measure the average total width for gamma decay of CN 
resonances located below neutron threshold. 

One of us (HAW) is grateful to P. Thirolf for drawing 
his attention to the problem, and for discussions. We are 
grateful to H. L. Harney, T. Papenbrock, and A. Richter 
for a reading of the manuscript and helpful comments. 
This work was partly supported through SFB634 by the 
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